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ABSTRACT 

A  beamspace  adaptive  beamformer  implementation 
for  the  rejection  of  cable  strum  self-noise  on  passive 
sonar  towed  arrays  is  presented.  The  approach  focuses 
on  the  implementation  of  a  white  noise  gain  constraint 
based  on  the  scaled  projection  technique  due  to  Cox  et 
al  [IEEE  Trans,  on  ASSP,  Vol.  35  (10),  Oct.  1987]. 
The  objective  is  to  balance  the  aggressive  adaptation 
necessary  for  nulling  the  strong  mainlobe  interference 
represented  by  cable  strum  against  the  conservative 
adaptation  required  for  protection  against  signal  self¬ 
nulling  associated  with  steering  vector  mismatch. 
Particular  attention  is  paid  to  the  definition  of  white 
noise  gain  as  the  metric  that  reflects  the  level  of 
mainlobe  adaptive  nulling  for  an  adaptive 
beamformer.  Adaptation  control  is  subsequently 
performed  through  the  implementation  of  a  constraint 
on  maximum  allowable  white  noise  gain  at  the  output 
of  the  adaptive  processor.  The  theoretical  development 
underlying  the  scaled  projection  based  constraint 
implementation  is  reviewed.  Towed  array  data  results 
depicting  the  performance  gain  of  the  new  ABF 
algorithm  optimized  for  strum  cancellation  relative  to 
that  of  a  more  conservative  baseline  ABF  algorithm 
are  presented. 

1.  INTRODUCTION 

Hydrodynamic  self-noise  on  passive  sonar  towed  arrays 
has  been  a  well-known  performance- limiting  factor  for 
ocean  acoustic  source  detection  at  low  frequency  [1]. 
High  wavenumber  mechanical  vibrations  are  induced  in 
the  array  by  vortex  shedding  associated  with 
hydrodynamic  flow  over  the  array  body  and  cable  scope. 
These  vibrations  are  know  to  couple  into  the  hydrophone 
array  as  coherent  acoustic  noise  sources  and  can  impair 
acoustic  detection  performance,  particularly  in  the 
forward  endfire  direction.  As  a  direct  consequence  of  its 
spatially  coherent  nature,  it  has  been  shown  that  cable 
stmm  noise  effects  can  be  mitigated  via  adaptive 


processing  [2].  In  this  work,  a  new  approach  to  coherent 
stmm  noise  mitigation,  based  on  a  beamspace  adaptive 
beamformer  (ABF)  architecture  with  a  white  noise  gain 
constraint  (WNGC)  that  emphasizes  mainlobe 
interference  nulling  is  introduced.  Finally,  data  results 
illustrating  the  performance  improvement  over  an  existing 
beamspace  ABF  algorithm  that  emphasizes  robustness  to 
mismatch-induced  self-nulling  are  presented. 

2.  THE  PHYSICS  OF  CABLE  STRUM 
2.1  Vortex  shedding 

When  an  array  is  subject  to  hydrodynamic  flow  with  a 
component  normal  to  its  axis,  a  wake  is  formed.  When  the 
velocity  of  the  transverse  flow  increases  beyond  a  certain 
threshold,  eddies,  or  vortices,  begin  to  form  and  separate 
from  the  wake.  Eventually  these  vortices  shed  from  the 
wake  in  an  asymmetric  fashion  [3].  This  asymmetric 
shedding  imparts  an  oscillatory  lift  force  locally  on  the 
array  which,  depending  on  the  properties  of  the  array  such 
as  tension  and  density,  can  excite  transverse  vibrations 
which  propagate  along  the  array  axis.  The  frequency  of 
vortex  shedding  in  hydrodynamic  flow  is  related  to 
properties  of  the  flow  and  the  array  via  the  empirically 
determined  Strouhal  relation  [1]: 


where  S  is  the  Strouhal  number,  equal  to  0.21  in  the 
laminar  flow  regime  characteristic  of  most  towed  array 
environments,  v  is  the  velocity  of  flow  normal  to  the  array 
axis,  and  d  is  the  cable  diameter.  Note  that  the  normal 
component  of  velocity  of  flow  can  vary  with  time  in 
response  to  platform  motion  and  local  inhomogeneities  in 
the  turbulent  medium. 

The  transfer  function  to  which  the  Strouhal  excitation  is 
applied  is  governed  by  the  wave  equation  subject  to  the 
boundary  conditions  of  the  array  under  tow.  For  example, 
assuming  fixed  boundary  conditions  for  the  array,  the 
preferred  frequencies  of  vibration  or  modes  of  the  array 
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corresponding  to  the  solution  of  the  wave  equation  is 
given  by: 


where  T  is  eable  tension,  nic  is  mass  per  unit  length  of  the  eable, 
and  L  is  the  eable  length.  Figure  1  depiets  notionally  the 
interaetion  of  the  Strouhal  exeitation  with  the  struetural  modes 
of  the  array.  Cable  strum  due  to  vortex  shedding  is  strongly 
exeited  when  the  Strouhal  exeitation  frequeney  is  elosely  aligned 
with  a  resonant  mode  of  the  eable  transfer  fimetion. 

2.2  Wavenumber-frequency  analysis 

The  deeomposition  of  an  array  snapshot  into  its  eonstituent 
aeoustie  and  non-aeoustie  eomponents  is  aeeomplished  using  a 
wavenumber-frequeney,  or  k-co,  transform.  The  k-co  transform  is 
a  2-d  FFT  in  spaee  and  time.  Maximum  unambiguous 
wavenumber  resolvable  is  equal  to  n/d,  where  d  is  the  sensor 
spaeing.  Resolution  in  wavenumber  is  governed  by  the  aperture 
length,  L.  For  non-dispersive  propagation,  frequeney  and 
wavenumber  are  linearly  related  via 


where  ep  equals  the  phase  speed  of  the  wavefront. 

Figures  2  and  3  depiet  k-co  plots  for  two  towed  arrays 
under  eonsideration  in  this  work.  The  first  exhibits 
superior  vibration  isolation  and  higher  resolution  due  to 
its  longer  aperture.  This  array  experienees  only  weak 
sidelobe  leakage  of  vibrational  modes  into  the  aeoustie 
eone.  As  sueh,  under  nominal  operating  eonditions,  this 
array  does  not  exhibit  a  pronouneed  eable  strum 
interferenee  problem.  The  second  array  is  characterized 
by  limited  vibration  isolation.  It  is  subjeet  to  signifieant 
leakage  of  vibrational  energy  into  the  aeoustie  eone  via 
mainlobe  penetration  in  forward  endfire.  Leakage  of 
vibrational  energy  into  aeoustie  forward  endfire  is  a 
strong  funetion  of  own- ship  tow  speed.  For  this  array, 
which  is  the  subjeet  array  for  this  paper,  eable  strum 
represents  a  signifieant  mainlobe  interferenee  problem. 


3.  BEAMS? ACE  ABF  FOR  CABLE 
STRUM 


The  beamspace  ABF  derives  its  eable  strum  nulling  eapability 
from  the  faet  that  near  endfire  the  beam  fan  is  partially 
eomposed  of  beams  steered  to  high  wavenumber  non-aeoustie 
spaee. 

For  eaeh  time  epoeh,  the  element  timeseries  are  transformed  to 
the  frequeney  domain  via  FFT.  A  beamspaee  eovarianee  matrix 
is  formed  for  eaeh  frequeney  bin  independently  and  a  7- 
dimensional  beamspaee  MVDR  weight  veetor  is  subsequently 
eomputed.  Adaptation  eontrol  is  governed  by  setting  a  limit  on 
the  maximum  allowable  white  noise  gain  for  the  adaptive 
proeessor. 

3.1  White  Noise  Gain 

White  noise  gain  (WNG)  is  defined  as  the  gain  applied  by 
the  adaptive  beamformer  to  a  spatially  white  input  noise 
proeess,  and  is  represented  by 

WNG  =  >v^>v, 

where  W  represents  the  MVDR  beamformer  steering 
veetor  given  by 

R~^v 

w  = - . 

v^R-^v 

The  veetor  v  represents  the  CBF  weight  veetor  and  the 
matrix  R  denotes  the  sample  eovarianee  for  the  eurrent 
proeessing  bin.  (Aetually,  the  beamformer  WNG  is  a 
quantity  equally  applieable  to  the  output  of  the  CBF 
beamformer,  expressed  as  V®V).  Beamformer  WNG  is  a 
measure  of  the  level  of  mainlobe  adaptive  nulling  effeeted 
by  the  beamformer  steering  veetor.  As  sueh,  a  eonstraint 
on  maximum  allowable  WNG  ean  be  used  to  eontrol  the 
level  of  mainlobe  adaptation  of  the  adaptive  beamformer 
relative  to  that  of  the  ideal  eonventional  beamformer: 

N 

Here  is  a  eonstant  ranging  from  1  to  infinity,  with  1 
representing  CBF  performanee  (no  adaptive  nulling 
eapability  and  best  robustness  to  mismateh)  and  infinity 
representing  MVDR  performanee  (most  adaptive  nulling 
eapability  and  most  sensitivity  to  mismateh).  Note  that 
under  this  eonvention,  the  quantity  1/N  represents  the 
WNG  of  the  eonventional  CBF  beamformer,  where  N 
equals  the  number  of  elements  in  the  array. 


The  ABF  arehiteeture  under  eonsideration  in  this  paper 
eonsists  of  a  frequeney-domain  beamspaee  adaptive 
beamformer.  The  adaptive  beamspaee  eonsists  of  a  7- 
dimensional  beam  fan  with  fixed  eosine  spaeing.  The 
beam  fan  translates  with  steering  direetion. 


The  relative  WNG  is  a  partieularly  important  metrie  to 
eonsider  when  the  souree  of  interferenee  lies  within  the 
beamformer  mainbeam.  Figure  4  depiets  the  behavior  of 
the  WNG  of  the  minimum  varianee  distortionless 


response  (MVDR)  ABF  relative  to  that  of  CBF  for  a 
simulation  seenario  in  whieh  an  interferer  is  swept  aeross 
eosine  spaee  and  permitted  to  penetrate  the  beamformer 
mainbeam.  An  elevation  in  WNG  results  from  the  ABF 
algorithm  attempting  to  drive  a  mainlobe  null  eoneurrent 
with  satisfying  the  MVDR  unity  gain  eonstraint  in  the 
steering  direetion.  The  three  inset  figures  show  ABF 
(shown  in  red)  and  CBF  (shown  in  blue)  beampattems  in 
the  vieinity  of  the  steering  direetion  for  three  different 
interferer  eosine  positions.  The  sequenee  attempts  to 
eonneet  the  WNG  eosine  dependenee  with  the  ABF 
beampattem  shape  as  the  interferer  eosine  approaehes  the 
steering  direetion.  When  the  interferer  is  far  in  the 
sidelobe  of  the  array  beampattem  (inset  3),  the  ABF  and 
CBF  mainlobe  beampattems  effeetively  overlap.  In  this 
ease,  a  simple  sidelobe  null  (not  pietured)  is  all  that  is 
needed  in  order  to  maximize  signal-to-interferenee-plus- 
noise  ratio  (SINR).  As  the  interferer  penetrates  the 
mainbeam,  a  squinting  or  splitting  of  the  adaptive 
beampattem  oeeurs  eoineident  with  the  introduetion  of  a 
mainlobe  null.  This  squinting  is  the  result  of  the 
beamformer ’s  attempt  to  maximize  SINR  by  trading  off 
interferenee  suppression  against  exeess  white  noise  gain 
in  the  vieinity  of  the  steering  direetion. 

3.2  Adaptivity/Robustness  Tradeoff 

Figure  4  illustrated  how  an  elevation  in  WNG  oeeurs  in 
response  to  a  mainlobe  interferer.  We  may  eonelude  that 
WNG  is  a  measure  of  the  mainbeam  adaptive  nulling 
being  performed  by  the  ABF.  It  is  important  to  understand 
that  the  ABF  algorithm  is  unable  to  distinguish  between 
most  forms  of  signal  model  mismateh  and  a  mainlobe 
interferer.  Thus,  the  ABF  will  interpret  steering  veetor 
mismateh  as  mainlobe  interferenee  and  attempt  to  eaneel 
it  as  well.  Some  degree  of  steering  veetor  mismateh  is 
unavoidable  in  real  towed  array  data  applieations. 
Common  sourees  of  mismateh  inelude  manifold 
uneertainty,  sensor  ealibration  error,  and  unmodeled 
multipath  propagation.  The  beamformer  signal  model  is 
based  on  an  assumption  of  a  perfeet  plane  wave  with 
known  sensor  gain  and  known  relative  sensor  loeation.  As 
the  ABF  algorithm  will  attempt  to  null  any  data 
eomponent  that  deviates  from  these  assumptions,  self¬ 
nulling  due  to  steering  veetor  mismateh  is  a  major 
eoneem.  By  imposing  a  eonstraint  on  the  maximum 
allowable  WNG  of  the  adaptive  beamformer,  robustness 
to  mismateh  indueed  nulling  may  be  introdueed. 

Analyses  of  towed  array  data  have  shown  that  to  effeet  a  useful 
level  of  strum  rejeetion  using  the  beamspaee  ABF  algorithm,  a 
fairly  aggressive  adaptation  strategy  is  required.  By  eontrast, 
signal  proteetion  against  self-nulling  in  the  eable  strum  band 
requires  a  very  eonservative  adaptation  approaeh.  In  this  work,  it 
was  empirieally  determined  that  a  WNGC  of  6  dB,  or  a 


maximum  allowable  WNG  of  4x  that  of  the  CBF  beamformer, 
represents  the  best  eompromise  between  mainlobe  eable  strum 
nulling  and  signal  preservation  in  the  presenee  of  mismateh. 


3.3  Adaptive  Weight  Power  Scaling 


The  white  noise  gain  eonstraint  (WNGC)  employed  in  the 
beamspaee  ABF  arehiteeture  is  based  on  the  sealed 
projeetion  teehnique  first  proposed  by  Cox  et  al  [4]. 


The  sealed  projeetion  WNGC  implementation  is 
eomposed  of  two  essential  parts.  First,  the  MVDR  weight 
veetor  is  deeomposed  into  two  orthogonal  eomponents, 
non-adaptive  and  adaptive  eomponents  respeetively, 
using  the  following  beamspaee  projeetion  operators: 


vv^ 

P  =  ———  P  =  T  —  P 

^  na  ^  ^  na 

V  V 


Seeond,  upon  a  WNG  threshold  exeeedanee,  the  adaptive 
eomponent  thus  isolated  is  sealed  sueh  that  the  WNGC  at 
the  beamformer  output  is  met  exaetly. 


The  orthogonal  deeomposition  prior  to  adaptive  weight 
sealing  is  important.  This  step  guarantees  that  the  weight 
sealing  will  be  applied  only  to  the  adaptive,  or  data- 
dependent,  eomponent  of  the  ABF  weight  veetor.  This 
insures  that  the  sealing  proeess  does  not  modify,  seale,  or 
rotate  the  beamformer  response  to  a  signal  that  is 
perfeetly  matehed  to  the  steering  veetor.  Consequently, 
the  sealing  preserves  the  eonstraint  of  distortionless 
response  in  the  steering  direetion.  The  adaptive 
eomponent  of  the  MVDR  weight  veetor  is  given  by: 

iF  =  P  w. 


It  is  straightforward  to  verify  that  the  non-adaptive  and 
adaptive  eomponents  derived  in  this  way  are  indeed 
orthogonal.  The  sealed  output  weight  veetor  is  then  given 
by: 


w  =  w  +  kw  , 

o  na 


where  the  sealar,  k,  represents  the  sealing  eoeffieient.  We 
then  speeifiy  the  WNGC  at  the  output  of  the  beamspaee 
ABF  proeessor  in  terms  of  a  multiplier  on  the  non- 
adaptive  WNG, 


w^T^Tw  <a  w^T^Tw  . 

rv  o  ^  —KA.  rv  ^  rr  . 


Here,  T  represents  the  7-dimensional  transformation  from 
element  spaee  to  adaptive  beamspaee.  For  a  6  dB  WNGC 
the  multiplier,  a,  is  equal  to  4.  Solving  the  eonstraint 
equation  results  in  a  quadratie  on  the  sealing  eoeffieient,  k 
[5].  The  result  is  two  solutions  for  k  whieh  meet  the 
eonstraint  exaetly.  We  ehoose  the  value  whieh  minimizes 
the  output  power  of  the  ABF.  This  proeedure  is  earried 
out  at  eaeh  proeessing  epoeh  and  for  eaeh  frequeney  bin 


independently.  A  geometrie  interpretation  of  the  weight 
sealing  proeedure  is  shown  in  Figure  5. 

4.  TOWED  ARRAY  DATA  RESULTS 

Figure  6  depicts  frequency-azimuth  (FRAZ)  plots  for  a 
typical  time  epoch  for  each  of  four  different  processors:  a) 
the  CBF  beamformer,  2)  the  conservative  baseline  ABF 
beamformer,  3)  the  aggressive  6  dB  WNGC  ABF 
optimized  for  strum  rejection,  and  4)  unconstrained 
MVDR.  The  baseline  ABF  represents  the  WNGC  as 
implemented  in  the  present  towed  array  processing 
system.  While  the  details  of  this  WNGC  implementation 
are  not  presented  here,  the  basic  design  philosophy  of  this 
ABF  algorithm  is  to  emphasize  robustness  to  mismatch 
effects.  Upon  exceeding  the  WNG  threshold,  set  in  the 
vicinity  of  2  dB,  the  baseline  ABF  scales  the  adaptive 
weight  vector  back  to  the  non-adaptive  or  CBF  weight 
vector.  This  severely  constrains  the  ability  of  the  baseline 
ABF  to  effectively  null  any  strong  mainlobe  interference 
such  as  cable  strum. 

In  Figure  6,  the  presenee  of  a  strong  interferenee  souree  with 
multiple  sidelobes  is  observed  near  broadside  in  the  CBF  FRAZ 
display.  As  expeeted  all  of  the  ABF  approaehes,  eonservative 
and  aggressive  alike,  demonstrate  the  eapaeity  to  null  sueh  a 
strong  diserete  sidelobe  interferenee  souree.  This  result  thus 
serves  as  a  useful  eonsisteney  eheek  of  algorithm 
implementation. 

Next,  we  direct  our  attention  to  the  cable  strum 
interference  near  forward  sector,  i.e.  near  cosine  equal  to 
1.  In  the  normalized  frequency  band  f=  0-0.3,  cable 
strum  is  observed  to  extend  over  a  wide  sector  of  cosine 
space  from  forward  endfire  to  near  broadside.  The 
important  differences  between  the  conservative  and 
aggressive  ABF  approaches  are  apparent  from  the  cable 
strum  rejection  performance  in  this  frequency  band.  The 
conservative  ABF  algorithm  does  very  little  to  reduce  the 
amplitude  of  the  strum  interference  in  forward  endfire. 
The  bearing  extent  of  the  strum  is  reduced  slightly.  With 
its  6  dB  WNGC,  the  bearing  extent  and  amplitude  of  the 
cable  strum  is  significantly  curtailed  relative  to  that  of  the 
conservative  baseline  ABF  algorithm. 

Figure  7  shows  raw  power  spectrum  density  plots  to 
further  illustrate  the  performance  improvement  realized 
with  increasingly  aggressive  adaptation.  Notice  that  the 
cable  strum  ABF  achieves  as  much  as  a  15  dB  local 
suppression  of  the  strum- dominated  noise  floor  in  the 
normalized  frequency  band  f  =  0.1 -0.3.  The  resulting 
noise  floor  suppression  uncovers  the  presence  of  a 
narrowband  feature  at  /  =  0.2  that  was  otherwise 


undetectable  in  the  CBF  and  baseline  ABF  configurations. 

Figure  8  shows  the  measured  WNG  plots  corresponding 
to  the  power  spectrum  density  plots  of  Figure  7.  The 
measured  WNG  illustrates  the  relationship  between 
WNGC  and  strum  rejection.  It  is  clear  that  at  a  WNGC  of 
6  dB  most  of  the  strum  noise  floor  suppression 
performance  is  realized.  Recall  that  the  point  here  is  to 
allow  the  ABF  algorithm  to  adapt  only  as  much  as 
necessary  to  effect  useful  cable  strum  noise  suppression. 

5.  CONCLUSIONS 

Mechanically  induced  towed  array  self-noise  limits 
detection  performance  in  passive  sonar  systems, 
particularly  at  forward  endfire.  In  this  work,  a  beamspace 
adaptive  beamforming  architecture  for  the  rejection  of 
strong  mainlobe  cable  strum  rejection  in  forward  endfire 
was  presented.  The  approach  focused  on  the  choice  of  a 
white  noise  gain  constraint  which  achieved  a  suitable 
balance  between  aggressive  adapation  for  effective  strum 
nulling  and  conservative  adaptation  for  robustness  to 
mismatch- induced  self-nulling.  A  WNGC  of  6  dB  relative 
to  the  WNG  for  the  non-adaptive  steering  vector  was 
empirically  determined  to  offer  the  best  balance.  The 
WNGC  implementation  was  based  on  the  scaled 
projection  technique  first  presented  by  Cox  et  al.  [4]. 
Significant  cable  strum  suppression  performance  was 
shown  to  be  possible,  on  the  order  of  15  dB  locally  within 
the  strum  interference  band. 
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Figure  1:  Notional  depiction  of  time  varying  Strouhal  excitation  and  array 
vibrational  modes. 
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Figure  2:  Frequency -wavenumber  plot  for  TB-29  array. 
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Figure  4:  Simulated  white  noise  gain  plot  and  associated  CBF 
and  MVDR  beampattems. 


Figure  5:  Geometric  interpretation  of  scaled  projection  weight 
scaling  technique. 


Figure  3:  Frequency -wavenumber  plot  for  TB-23  array. 
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Figure  6:  Frequency -Azimuth 
plots  depicting  reduction  in 
cable  strum  cosine  extent  with 
increasingly  more  aggressive 
adaptation. 
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Figure  7:  Raw  power 
spectrum  density  plots, 
illustrating  strum  noise  floor 
suppression  with  increasing 
levels  of  adaptation.  Locally 
as  much  as  15  dB 
improvement  in  strum  noise 
suppression  is  realized  with 
6  dB  WNGC. 
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Figure  8:  Measured  WNG 
plots  corresponding  to  time 
epoch  that  produced  PSD 
plots  in  Figure  7.  Most  of 
the  performance  benefit  is 
reahzed  with  the  WNGC 
eaual  to  6  dB. 
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